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ABSTRACT 


The  results  of  a  previous  study  of  the  concept  of  a  stirred  reactor 
laser  with  lasing  restricted  to  the  zeroth  fundamental  band  are  compared 
with  the  predicted  power  output  of  a  multilevel  laser  for  the  hydrogen- 
fluorine- helium  (H2_F2"He)  system.  Lasing  is  assumed  to  occur  on 
all  fundamental  bands  for  which  the  calculated  maximum  gain  equals 
the  cavity  threshold  value.  Power  output  is  significantly  increased 
over  the  zeroth  band  calculations,  and  the  effect  of  cavity  pressure  is 
shown  to  be  less  pronounced  because  of  the  cascade  enhancement  from 
higher  bands.  Vibrational  de-excitation  caused  by  collisional  processes 
appears  to  limit  this  kind  of  chemical  system  to  low  cavity  pressures. 
The  effect  of  cavity  losses  was  evaluated  with  three  assumed  values 
for  the  cavity  loss  coefficient.  These  values  were  estimated  as  a 
small  fraction  (0.  1)  of  the  average  small  signal  gain  coefficient  at  dif¬ 
ferent  cavity  residence  times. 
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NOMENCLATURE 

Einstein  spontaneous  emission  coefficient 
Mirror  absorption  loss  coefficient 
Vacuum  speed  of  light 

Constant  in  rotational  energy  expression 

Gain  coefficient 

Small  signal- gain  coefficient 

Degeneracy 

Planck's  constant 

Intensity 

Rotational  quantum  number 

Degree  of  inhomogeneous  broadening 

Boltzmann's  constant 

Length  of  active  medium  in  cavity 

Mass  flow  rate  into  reactor 
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Pv 


Number  of  vibrational  levels 
Refractive  index 
Cavity  pressure 
Power  output 
Useful  power  output 

Rotational  partition  function  for  vibrational  level  v 

Stimulated  emission  function 

Mirror  reflection  coefficient 

Integrated  line  intensity 

Translational  temperature 

Cavity  residence  time 

Volume  of  reactor 

Vibrational  quantum  number 

Molecular  weight 

Energy 

Saturation  parameter 
Distance  along  cavity 
Coefficients  in  analytical  model 
Loss  co efficient 

Mole/ mass  ratio  of  vibration  level  v 
Chemical  production  rate  of  vibrational  level  v 

Stimulated  emission/absorption  production  rate  of 
vibrational  level  v 

Freauency 

Frequency  at  line  center 

Local  value  of  gain 

Gas  density 

Photon  flux  density 

Stimulated  emission  rate  constant 
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rv 

Gain 

u 

Wave  number 

SUPERSCRIPTS 

D 

Denotes  Doppler  effect 

o 

Denotes  reactor  feed  quantity 

V 

Denotes  Voigt  effect 

± 

Direction  indicator 

— 

Average  quantity 

SUBSCRIPTS 

D 

Denotes  Doppler  effect 

opt 

Optimum 

v,  o,  1,  etc. 

Denotes  vibrational  levels 

V 

Denotes  frequency 
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SECTION  I 
INTRODUCTION 

The  concept  of  a  continuous  wave  (cw)  stirred  reactor  chemical 
laser  was  presented  in  Ref.  1  for  a  hydrogen/fluorine/helium  (H2-F2- 
He)  system.  Calculations  were  carried  out  with  stimulated  emission/ 
absorption  transitions  restricted  to  the  zeroth  fundamental  vibrational 
band.  An  optical  gain  coefficient  was  calculated  for  the  rotational  line 
having  maximum  gain,  and  lasing  was  assumed  to  occur  at  the  center 
of  the  Doppler-broadened  line.  The  single  fundamental  band  restric¬ 
tion  was  used  so  that  one  could  specify  a  radiation  flux  density  for 
parametric  calculations  without  having  to  specify  detailed  optical  cav¬ 
ity  boundary  conditions  and  losses.  The  results  correspond  to  a  laser 
cavity  with  a  mirror  at  one  end  and  a  grating  at  the  other  to  select  the 
zeroth  fundamental  band.  In  this  report,  all  fundamental  bands  (up  to 
six)  with  sufficient  optical  gain  were  allowed  to  lase.  A  cavity  gain 
threshold  equation  was  used  to  determine  the  radiation  flux  densities 
and  which  bands  were  capable  of  lasing.  The  predicted  power  output 
reported  in  Ref.  1  indicated  a  dramatic  decrease  as  the  cavity  pres¬ 
sure  level  was  increased.  As  cavity  pressure  increases,  collisional 
relaxation  rates  increase,  thus  reducing  the  optical  gain  coefficient. 
Allowing  all  fundamental  bands  to  lase  should  reduce  this  effect  since 
the  decrease  in  gain  on  higher  bands  is  offset  by  an  equivalent  increase 
in  gain  on  lower  bands.  This  cascade  enhancement  is  particularly  preva¬ 
lent  during  lasing  where  the  reduction  in  population  of  an  upper  vibra¬ 
tional  level  appears  as  an  increase  in  population  of  the  next  lower  level. 

The  purposes  of  this  research  were  (1)  to  apply  the  AEDC  expertise 
in  propulsion  technology  and  reaction  kinetics  to  the  stirred  reactor 
concept  in  order  to  determine  its  potential  as  a  chemical  laser  device 
and  (2)  to  become  familiar  with  chemical  laser  devices  which  could 
conceivably  be  tested  in  AEDC  facilities. 

A  simplified  analysis  of  the  effect  of  cavity  losses  is  included  to 
determine  the  effect  on  power  output  and  to  determine  whether  an  opti¬ 
mum  mirror  transmission  coefficient  exists  as  predicted  for  other 
chemical  laser  models  (Ref.  2). 

SECTION  II 

THEORETICAL  MODEL 

2.1  CHEMICAL  AND  RADIATIVE  INTERCHANGE  EQUATIONS 

in  accordance  with  Ref.  1,  the  vibration  levels  v  =  0,  1,  --n  of 
the  radiating  molecule  are  considered  as  an  individual  chemical 
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species  denoted  by  yv.  The  equations  which  describe  the  chemical  and 
radiative  interchanges  (assuming  only  fundamental  band  transitions  are 
allowed)  in  a  perfectly  stirred  reactor  are 


-  yv  =  -tcyv  —  and  v  =  °’n  (2.  1) 

C  8 

O 

where  yv  is  the  concentration  (moles/gm)  in  the  reactor  feed,  yv  is  the 
concentration  in  the  reactor,  7Vc  is  the  rate  of  production  of  vibrational 

level  v  from  chemical  reaction,  7V  is  the  rate  of  production  attribut- 

s 

able  to  stimulated  emission/absorption  transitions  (spontaneous  radia¬ 
tion  assumed  to  be  negligible),  and  tc  is  the  average  residence  time  of 
particles  in  the  reactor.  The  quantity  7Vc  is  evaluated  conventionally 

in  terms  of  a  forward- rate  constant  and  a  backward- rate  constant 

(obtained  from  the  forward- rate  constant  and  an  equilibrium  constant). 

The  quantity  yv  is  evaluated  as  the  net  production  rate  of  vibrational 
s 

level  v  as  follows: 


Vv  ~~Pv  v  “  1v—  iV v-i^  Pv  (yv+i  “  IvlO 

s  v— 1  v 

Pvy  =  ^j^~3  ^^vUv-De.pf-Fv^Uv-D^/QS^j 


(2.2) 

(2.3) 


qv  =  «P  {[Fv+iUv-D  -  Fv(Jv)]  (2.  4) 

The  quantity  pv  is  the  radiation  flux  density,  qv  is  the  minimum 
population  ratio  for  lasing  to  occur,  is  the  Einstein  coefficient 

for  spontaneous  transitions,  gv(Jv  -  1)  is  the  degeneracy  of  the  tran¬ 
sition,  Jv  denotes  the  rotational  line  having  maximum  gain,  Fv(Jv)hc/k 
is  a  characteristic  rotational  temperature,  is  the  rotational  parti¬ 
tion  function  for  vibrational  level,  v,  wv  is  the  wave  number  of  the 
transition,  h  is  Planck's  constant,  k  is  Boltzmann's  constant,  c  is  the 
vacuum  speed  of  light,  and  T  is  the  rotational  temperature  (assumed 
equal  to  the  translational  temperature). 


2.2  CAVITY  THRESHOLD  EQUATION  AND  RADIATION  FLUX  DENSITIES 

The  problem  in  solving  Eqs.  (2.  1)  is  that  of  specifying  the  radia¬ 
tion  flux  densities,  pv^.  In  Ref.  1,  only  pv  was  considered,  with 
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radiation  transitions  between  vibrational  levels  v  =  1  and  v  =  0  allowed. 
Values  of  pv  were  assumed  and  parametric  calculations  carried  out  for 
various  pressure  levels  and  residence  times,  tc.  The  radiation  flux 
densities  may  be  calculated  implicitly  through  some  algebraic  manipu¬ 
lation  of  Eqs.  (2.  1)  and  the  introduction  of  the  laser  (gain)  threshold 
equation  (Ref.  2). 


2Lgv  =  2L/9l  -  £n  r,r0  (2.5) 

where  L  is  the  active  cavity  length,  gv  is  the  gain  coefficient,  is  a 
cavity  loss  coefficient,  and  r^  and  r2  are  the  reflection  coefficients  of 
the  cavity  mirrors. 

Substituting  Eq.  (2.  2)  in  Eqs.  (2.  1)  and  noting  that  ~pv  .  is  zero, 
one  obtains 


y 


o 


(2.  6) 


and 


y“  -  y» 


n 


2 

v=0 


yv 


(2.  7) 


There  is  thus  one  equation  instead  of  n,  and  an  additional  n-1  equa¬ 
tions  are  required  for  the  remaining  n-1  unknown  yv.  Equation  (2.  5) 
together  with  the  spectral  line  properties  for  each  transition  is  used  to 
provide  the  additional  equations. 

The  gain  coefficient  gv  is  related  to  the  spectral  line  through  the 
integrated  line  intensity  Sv  as 


The  coefficient  is  evaluated  in  Ref.  1  as  a  Voigt  broadening 
coefficient  in  order  to  account  for  combined  collision/Doppler  broaden¬ 
ing  processes,  gP  is  the  pure  Doppler-broadened  gain  coefficient,  and 
Ai/j}  is  the  Doppler  half- width  at  half  maximum  intensity  for  the  transi¬ 
tion  at  frequency  v. 
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From  the  definition  of  Sy  it  can  be  determined  that 

-Sv  =  htuvPNAaC(>'v+l- 


where 


8ffhai i 


•C1  gv(Jv-l)  exp 


From  Eqs.  (2.  5),  {2.  8),  and  (2.  9)  one  can  derive 


yv+i  =  +  av  and  v=  °*n-1 


where 


KVhN . 


[mT  J_ 

1  *  P  < 


(2.  9) 


(2.  10) 


(2.  11) 


(2.  12) 


Here,  is  Avagadro's  number,  p  is  the  mixture  density,  and  W  is  the 
molecular  weight  of  the  radiating  molecule. 

Equations  (2.  11)  thus  enable  yv,  v  =  1,  n  to  be  evaluated  in  terms 
of  t0-  Nonradiating  species  are  assigned  an  equation  similar  to  Eq. 

(2.  1)  (without  the  radiation  production  team),  except  for  those  evaluated 
with  atom  conservation  equations.  The  system  of  equations  is  solved 
by  the  iterative  technique  described  in  Refs.  1  and  3  after  a  loss  coef¬ 
ficient  is  defined.  The  mixture  density  and  temperature  are  obtained 
from  the  equation  of  state  and  an  energy  balance  which  includes  the 
radiated  power  and  radiative  losses. 


SECTION  III 

LASER  POWER  OUTPUT  WITH  CAVITY  LOSSES 


The  intensity  in  a  laser  cavity  with  mirrors  at  z  =  0  and  z  =  L  is 
related  to  the  gain  coefficient  as  follows  (losses  are  included  later): 

dl+  ar 

aT  “  ^  and  TT  "  <3*  D 

In  a  perfectly  stirred  reactor,  the  gain  coefficient  gv  is  independ¬ 
ent  of  position  z;  thus,  Eq.  (3.  1)  can  be  integrated  to  give 
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I+(z)  =  efl+(l)  and  £(z)  =  i;(2)e^"^  (3.  2) 

where 


and  rv  =  g^L 

Under  lasing  conditions,  the  intensity  has  a  very  narrow  spike  at  vQ; 
thus,  Eqs.  (3.  2)  can  be  integrated  over  the  narrow  frequency  band¬ 
width  for  the  spike  to  give 

£  r  —  £ 

I+(z)  =  e  ‘'“I+U)  and  I“(2)  =  e  "°  Vo  H2)  (3.3) 

At  the  mirror  boundaries, 

I+(l)  =  rjHl)  and  I“(2)  -=  ^^(2)  (3.  4) 


1 


rlr2 


(3.  5) 


r+d)  = 
r(2) 


(3.6) 


Equatiop  (3.  5)  is  the  cavity  threshold  equation  describing  the  radiative 
steady- state  condition  that  says  gain  equals  loss.  In  other  words,  the 
increase  in  intensity  for  two  passes  through  the  medium  is  balanced  by 
the  output  or  losses  at  the  two  mirrors. 


The  intensities  at  the  boundaries  can  be  determined  by  calculating 
an  average  intensity  I.  Thus, 


V-  =  I  ;L  rtzMz  =  J_  /  ^  l%  )d$, 

n  r*j  O  o 


(3.  7) 


I+(l),  and  T 


I- (2) 


(3.  8) 
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The  average  intensity  is  given  by 

T  =  T+  +  r  =  e  1  [i+(D  +  r(2)] 


(3.  9) 


3.1  OUTPUT  POWER 

The  total  output  power  is  given  by  the  product  of  the  energy  per 
photon,  the  number  of  photons  per  unit  volume  to  leave  the  cavity,  and 
the  cavity  volume,  as  follows: 

PR  =  2  hcovPyAPl/^(yw-  q^y^V  (3.10) 

The  output  power  is  also  given  by  the  flux  crossing  the  mirror  surfaces, 
as  shown  by 


PR  =  [(1  -  rjH-Q)  +  (1  -  r2)l+(2)l  ^ 


(3.  11) 


It  may  be  assumed  that 


ri  ■  r2 


r 


m 


.■.pR  -  ♦  na]  l 

m 


2(1-  r  ) 
m 


r 

m 


2d  -  r  ) 
m 


r 

m 


o 


(3.  12) 

(3.  13) 
(3.  14) 

(3.  151 

(3.  16) 
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Equation  (3.  16)  relates  the  total  power  output  to  the  average  cavity 
intensity 


3.2  CAVITY  POWER  OUTPUT  WITH  LOSSES 

It  is  necessary  to  add  j3j_I±  to  t^le  right-hand  side  of  the  radiative 
transfer  equations  (3.  1);  thus,  the  rv  are  replaced  by  tv  -  dj^L,  and 
Eq.  (3.  5)  becomes 


-  -6*  +  )3lL  (3.  17) 

The  useful  power  output  is  readily  evaluated  as 

(1-^)P"  <3'18) 

The  quantity  tv  is  the  total  output  of  which  dLL  is  lost  to  absorption, 
scattering,  and  diffraction  losses. 

Inspection  of  Eq.  {3.  18)  indicates  that  a  maximum  useful  power 
output  must  occur  at  some  appropriate  ratio  of  j3L/ gv  since  power  out¬ 
put  is  zero  if  I  *  0  and  also  if  gv  =  0l  (corresponding  to  r  *  1  in  Eq. 

(3.  17)).  Since  Pr/hi  is  a  function  of  j3L  and  gv,  it  is  not  possible  to 
evaluate  the  maximum  Pj^/m  analytically  from  Eq.  (3.  18).  Under 

certain  circumstances  an  analytical  solution  for  laser  power  output 
can  be  obtained  (Refs.  4  through  9).  The  analyses  (based  on  gain  satu¬ 
ration  formulas  given  in  Ref.  5)  predict  that  an  optimum  output  coupling 
exists  for  maximum  power  output  when  cavity  losses  are  included  in  the 
analysis.  This  is  in  contrast  to  the  no-loss  case  in  which  the  radiation 
flux  density  in  the  cavity  can  be  increased  to  infinity  as  shown  in  Ref.  1. 
Two  simple  formulas  can  be  derived  for  the  optimum  output  coupling 
for  the  following  limiting  cases: 

1.  Homogeneous  (pure  collision)  broadening 

-(fr)  -  '/eJl  -  ft.  (3. 19) 

N  'opt 

2.  Inhomogeneous  (pure  Doppler)  broadening 

(3-2o) 


7 


AE  DC -TR -73-22 


It  is  not  immediately  apparent  what  saturation  formula  applies  to 
the  chemical  laser  described  in  this  report.  Equations  (3.  19)  and  (3.20) 
are  limiting  cases  derived  from  an  analysis  of  a  helium/neon  (He-Ne) 
laser  (Ref.  10).  For  this  laser,  production  terms  corresponding  to 
the  right-hand  side  of  Eqs.  (2.  1)  are  steady- state  values  and.  hence, 
are  equal  to  zero.  Given  these  conditions,  one  can  find  a  simple  solution 
for  the  lasing  species  concentration  by  algebraic  manipulation. 


3.3  MIRROR  ABSORPTION  LOSSES 

The  previous  development  provided  a  means  for  estimating  the  in¬ 
ternal  cavity  losses.  Mirror  surface  losses  are  accounted  for  by  re¬ 
ducing  the  value  of  r  to  its  actual  value  and  including  a  mirror  absorp¬ 
tion  coefficient,  a.  The  actual  value  of  the  transmissivity  of  each  mir¬ 
ror  is  then  given  by  1-r-a. 

Equations  (3.4)  through  (3.6)  become 

I+(l)  =  rr(l)  and  I~(2)  =  rI+(2>  (3.21) 


1  =  re  0 

(3.22) 

and 

I+(l)=  r(2) 

(3.23) 

Equation  (3.  11)  becomes 

PR  =  [( 1 — r— a)I  (1)  +  (l-rwa)I+(2)]  1 

u 

(3.24) 

=  (1— r-a)2I“(l)  J 

(3.  25) 

=  (i-r-a)  2i+(D  y 

r  L 

(3.  26) 

(3.27.) 

Thus,  with  mirror  absorption  losses  an  extra  correction,  \ 

f'j-r-aX  . 

a=r)- 
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introduced.  This  correction  factor  is  very  useful  if  closed-cavity 
experiments  are  conducted.  If  1  =  r  +  a,  the  cavity  power  is  lost  into 
the  mirror,  which  may  be  water-cooled  for  measurement  of  the  flux. 

The  value  of  j3L  to  be  used  cannot  be  readily  evaluated  theoretically. 
In  a  perfectly  stirred  reactor  laser,  the  losses  would  be  primarily 
scattering  and  diffraction  losses  attributable  to  the  local  fluctuations 
generated  by  the  intense,  fine- scale  turbulence.  In  a  well- stirred  re¬ 
actor,  in  which  the  mixing  is  not  entirely  complete,  some  absorption 
would  occur  also,  since  all  the  eddies  would  not  have  undergone  exactly 
the  same  time  history  in  the  reactor.  Large  absorption  losses  would 
occur  if  the  eddies  were  large;  i.  e. ,  mixing  occurs  on  a  macroscale 
and  not  on  a  molecular  scale.  In  this  latter  case,  the  ensemble  of  ed¬ 
dies  is  considered  as  a  collection  of  individual  stirred  reactors  which 
remain  segregated.  The  residence  time  distribution  in  a  stirred  re¬ 
actor  is  given  by  exp(-t/tc)/tc  where  tc  is  a  mean  residence  time  in 
the  reactor.  The  output  concentration  in  the  segregated  reactor  is 
given  by 

yv  =  r-  f  yy(0  exp  (— t/tc)dt 

c  ° 

where  7v(t)  is  obtained  from  the  plug  flow  form  of  Eq.  (2.  1): 

dy  v 

■Jr  •  ♦  n. 


Segregation  of  fluid  elements  in  the  reactor  is  possible  if  there  is 
laminar  or  transitional  flow  in  the  reactor.  This  possibility  may  be 
very  real  at  cavity  pressures  much  below  0.  1  atm.  To  the  author's 
knowledge,  experimental  data  from  reactors  at  such  low  pressures 
has  not  been  obtained. 

An  analysis  of  the  effects  of  incomplete  molecular  scale  mixing  at 
atmospheric  pressure  is  given  in  Ref.  11  for  anthracite  flames  in  a 
swirl  combustor. 

Current  theoretical  models  of  (cw)  chemical  lasers  do  not  provide 
a  means  of  caluclating  Pj^;  thus,  a  simplified  method  of  accounting  for 
these  loss  mechanisms  is  used.  Obviously  the  value  of  Pj^  must  be  less 
than  the  small  signal  gain  and  for  this  reason,  values  of  /3l  =  0*  001, 

0.  0025,  and  0.  005  were  chosen  to  take  account  of  significant  absorption 
losses.  It  will  be  seen  in  the  next  section  that  these  values  correspond 
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to  one  tenth  of  the  average  small  signal  gain  coefficient  at  different 
values  of  residence  time,  tc. 


SECTION  IV 

RESULTS  AND  DISCUSSION 


Results  of  calculations  are  presented  for  a  H2-F2-He  system  in 
the  molar  ratios  2:1:10.  A  preburner  temperature  is  assumed  such 
that  the  temperature  of  the  fluorine  molecules  and  fluorine  atoms  in 
the  reactor  feed  is  600°K  at  the  cavity  pressure.  Thus,  as  the  cavity 
pressure  is  increased,  the  equilibrium  fluorine  atom  concentration 
decreases.  Hydrogen  and  helium  are  assumed  to  be  at  300°K  in  the 
reactor  feed.  The  fixed  preburner  temperature  is  considered  a  rea¬ 
sonable  limit  for  the  integrity  of  combustor  wall  materials.  Refrac¬ 
tive  index  fluctuations  in  the  subsonic  cavity  flow  are  considered  to  be 
negligible  as  discussed  in  Ref.  1.  Some  results  from  Ref.  1  are  in¬ 
cluded  to  illustrate  the  enhancement  caused  by  the  cascade  effect  of 
depopulating  all  fundamental  bands  capable  of  lasing.  Results  with  and 
without  cavity  losses  are  also  shown,  and  only  one  axial  mode  at  each 
line  center  is  considered. 

The  predicted  power  output  when  emission/ absorption  is  restricted 
tov=l^v=0is  shown  in  Figs.  1  through  3,  Appendix.  The  power 
output  is  plotted  as  a  function  of  mass  flow  rate  per  unit  volume,  m/V. 
The  average  residence  time  of  particles  in  the  reactor  is  simply 
tc  *  pV/m,  or,  since  p  is  approximately  constant  for  a  given  cavity 
pressure,  tc«  V/m.  The  maximum  power  output  is  readily  obtained 
by  setting  av  =  0  in  Eq.  (2.  12);  i.  e.,  the  radiation  flux  density  is  suf¬ 
ficiently  high  to  force  71  to  the  minimum  value  for  lasing:-  y^  =  q070 
(pv — ►«*>  as  go-—*-0).  The  power  output  for  =  0  with  all  fundamental 
bands  allowed  to  lase  (providing  the  threshold  equation  (2.  5)  is  satis¬ 
fied)  is  shown  in  Figs.  4  through  6.  The  cascade  enhancement  signifi¬ 
cantly  increases  the  output  over  single-band  calculations  shown  in  Figs. 
1  through  3.  The  reduction  in  peak  power  output  as  pressure  increases 
is  not  so  dramatic  in  Fig.  4  as  it  is  in  Fig.  2.  The  peak  is  readily  ex¬ 
plained  by  a  consideration  of  the  combined  actions  of  chemical  pumping, 
collisional  deactivation,  and  stimulated  radiation  depopulation  process¬ 
es.  At  very  short  times  tc  (large  rh/V),  chemical  pumping  is  attempt¬ 
ing  to  build  up  the  population  while  radiation  depopulation  reduces  it. 

As  tc  increases,  population  levels  increase  and  collisional  deactivation 
becomes  significant.  Collisional  deactivation  results  in  a  higher 
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mixture  temperature,  which  shifts  the  maximum  gain  line  to  higher  J 
values,  in  addition  to  the  reduction  in  population  ratios  of  adjacent  vib¬ 
rational  levels.  This  "loss"  is  in  contrast  to  the  carbon  dioxide  (CO2) 
gas-dynamic  laser  behavior  where  vibrational  levels  of  the  Vi  and  1^2 
modes  are  preferentially  reduced  by  rapid  collisional  deactivation  to 
maintain  adequate  gain  between  the  1/3  and  modes.  As  the  pressure 
is  increased,  collisional  deactivation  becomes  significant  earlier  be¬ 
cause  of  increased  population  levels,  and  the  peak  power  shifts  to 
shorter  and  shorter  times  {chemical  pumping  rates  also  increase  be¬ 
cause  of  increased  pressures).  The  very  short  residence  times  at 
pressures  greater  than  0.  025  atm  probably  preclude  higher  pressure 
operation  because  of  the  inability  to  mix  reactants  fast  enough  to  sat¬ 
isfy  the  reactor  mixing  criterion- -mixing  time  less  than  residence 
time. 

The  effect  of  introducing  cavity  losses  is  shown  in  Figs  7  and  8. 
The  power  output  is  significantly  reduced.  Values  of  the  loss  coeffi¬ 
cient,  were  chosen  from  the  small  signal  gain  calculations  shown  in 
Figs.  9  and  10  {corresponding  population  ratios  are  shown  in  Fig.  11). 
Values  of  j3j__  =  0.  1  gQ  with  gQ  =  0.  05,  0.  025,  and  0.  01  were  selected 
as  representative.  A  constant  average  value  of  gQ  was  used  for  con¬ 
venience  and  should  provide  a  conservative  estimate  of  the  useful 
power  output  of  a  well-stirred  reactor.  The  values  of  J3L  used  result 
in  a  significant  reduction  in  predicted  useful  power  output  with  an 
optimum  transmission  coefficient  r  «  0.  8.  The  results  shown  in  Fig. 

8  correspond  to  a  predominantly  Doppler- broadened  gain  coefficient; 
however,  the  value  of  r  =  0.  8  would  lead  to  a  negative  value  of  gD  if 
Eq.  {3.2)  were  applied,  since  |3l  <  -inr/2L;  hence,  this  equation  does 
not  apply  as  expected.  The  effect  of  mirror  absorption  losses  is 
shown  in  Fig.  12  for  a  =  0.  01  and  for  j3 l  ■  0. 

4.1  CAVITY  RESIDENCE  TIME  AND  MIXING 

The  cavity  residence  times  for  maximum  power  operation  are 
quite  short  {tc  *  100  psec  for  pressures  of  20  Torr  or  less);  however, 
if  there  is  adequate  turbulence  at  these  low  pressure  levels,  the  time 
required  to  mix  the  reactants  may  still  be  sufficiently  less  than  this  to 
satisfy  the  stirred  reactor  assumption.  Experiments  with  stirred  re¬ 
actors  have  been  carried  out  at  pressures  as  low  as  0.  1  atm  with  resi¬ 
dence  times  on  the  order  of  40  to  60  Msec  (Ref.  12). 
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Supersonic  flow  lasers  reported  in  the  literature  are  known  to  be 
diffusion  limited.  Supersonic  velocities  are  required  (1)  to  provide 
sufficiently  high  particle  transit  times  and  finite  beam  lengths  in  the 
flow  direction,  and  (2)  to  provide  some  diffuser  pressure  recovery  to 
reduce  exhaust  gas  pumping  requirements.  The  basic  stirred  reactor 
assumption  is  that  the  process  is  chemical- kinetic  limited.  The  mix¬ 
ing  time  must,  therefore,  be  short  compared  with  reactor  residence 
times.  It  has  been  well  established  (Refs.  13  and  14)  that  the  mixing 
rates  in  turbulent  shear  flows  are  determined  by  the  conversion  of 
mean  velocity  (U)  to  the  fluctuating  velocity  (u').  That  is,  the  mixing 

time,  tm,  is  proportional  to  the  ratio  of  these  velocities: 

/ 

U 

t  oe  - 

™  u  , 

a.  Incompressible  jet  into  still  air,  (— )  <  0.  2 

u  max 

b.  Incompressible  jet  into  coflowing  stream,  £=)  <  0.  2 

u  max 

c.  Supersonic  jet  into  coflowing  stream,  (jri  <  0.  05 

u  max 

u 9 

d.  Stirred  reactor,  (- )  >  1 

u  max 

The  above  numbers  have  been  established  in  experiments  at  pres¬ 
sures  of  0.  1  atmospheres  or  higher.  No  information  exists  in  the  lit¬ 
erature  for  pressures  lower  than  these.  The  pressure  levels  of  inter¬ 
est  for  chemical  lasers  are  significantly  less  than  this  -  on  the  order 
of  a  few  Torr  in  supersonic  flow  lasers.  Since  the  Reynolds  number 
is  inversely  proportional  to  pressure,  it  is  expected  that  (u7U)max 
will  be  even  lower  than  those  given  above.  In  fact,  it  is  possible  that 
the  supersonic  flow  lasers  have  transitional  or  laminar  mixing  rates. 

The  mixing  rate  in  a  stirred  reactor  is  usually  a  factor  of  20  to  50 
times  greater  than  that  for  the  coflowing,  supersonic,  diffusion  flame. 

In  addition,  the  time  constant  for  the  degree  of  chemical  conversion 
in  a  stirred  reactor  can  be  significantly  longer  than  that  in  a  premixed 
flame  (plug  flow  reactor)  as  shown  in  Ref.  3. 

For  these  reasons,  it  is  hoped  that  chemical-kinetically  controlled 
processes  in  a  stirred  reactor  laser  are  possible.  Only  experiments 
can  determine  if  there  is  adequate  turbulence  to  satisfy  the  basic  stirred 
reactor  assumption. 
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4.2  NUMERICAL  PROBLEMS 

Some  difficulties  were  experienced  in  obtaining  solutions  at  long 
residence  times  because  of  the  method  used  to  eliminate  lasing  be¬ 
tween  adjacent  lines  since  Eq.  (2.  5)  was  no  longer  satisfied.  The  cur¬ 
rent  version  of  the  program  assumes  lasing  is  terminated  on  succes¬ 
sively  higher  vibrational  levels  as  the  cascade  effect  fills  up  lower 
levels.  In  fact,  the  gain  is  not  distributed  sequentially  as  shown  in 
Figs.  8  and  9  even  for  no  power  output.  This  is,  of  course,  caused  by 
the  nonuniform  pumping  of  the  chemical  reactions  (Ref.  15). 


SECTION  V 
CONCLUSIONS 


A  method  of  predicting  the  power  output  from  vibrationally  excited 
molecules  in  a  stirred  reactor  laser  has  been  developed  and  the  follow¬ 
ing  conclusions  have  been  reached: 

1.  Radiation  flux  densities  are  determined  by 
the  requirement  that  the  gain  between  vi¬ 
brational  levels  satisfy  the  threshold  equa¬ 
tion  for  a  resonance  cavity. 

2.  The  power  output  for  multilevel  lasing  is 
significantly  greater  than  that  for  the  sim¬ 
ple  two- level  laser  analyzed  previously. 

3.  The  effect  of  cavity  pressure  level  is  less 
pronounced  for  the  multilevel  than  for  the 
two- level  laser;  however,  the  H2-F2-He 
system  appears  to  be  limited  to  low  pres¬ 
sure  levels  because  of  the  influence  of 
collisional  deactivation  processes. 

4.  Losses  reduce  the  power  output  signifi¬ 
cantly,  with  an  optimum  mirror  reflec¬ 
tion  coefficient  r  =  0.  8  for  the  chemical 
laser  analyzed  in  this  report. 

5.  Cavity  residence  times  are  quite  short 
(tc  =  100  /nsec  for  pressures  of  20  Torr 
or  less),  but  if  there  is  adequate  turbu¬ 
lence  for  mixing  on  a  molecular  scale. 
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the  mixing  time  may  still  be  sufficiently 
less  than  this  to  satisfy  the  stirred  re¬ 
actor  assumption. 
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APPENDIX 

ILLUSTRATIONS 
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Fig.  1  Specific  Power  versus  Mass  Flow  Rate  per  Unit  Volume:  Effect  of  Radiation  Flux  Density 
at  Constant  Pressure  and  Enthalpy 


Specific  Power,  — ,  w/g  sec 


V 

Fig.  2  Specific  Power  versus  Mass  Flow  Bate  per  Unit  Volume:  Effect  of  Pressure  at  Constant  Enthalpy 
and  Radiation  Flux  Density 
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Fig.  4  Specific  Power  versus  Mass  Flow  Rate  per  Unit  Volume:  Effect  of  Pressure  at  Constant  (Low)  Gain 
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a.  0L  =  0.001 

Fig.  7  Specific  Power  versus  Residence  Time:  Effect  of  Mirror  Reflection  Coefficient  at 
a  Cavity  Pressure  of  0.02  atm 
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c.  0L  =  0.005 
Fig.  7  Concluded 
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Useful  Specific  Power,  - ,  w/g  sec 
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a.  0L  =  0.001 

Fig.  8  Useful  Specific  Power  versus  Mirror  Reflection  Coefficient  at 
a  Cavity  Pressure  of  0.02  atm 
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Fig.  11  Inversion  Ratio  versus  Residence  Time  at  Zero  Power 
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Fig.  12  Effect  of  Mirror  Absorption  Losses  on  Useful  Specific  Power 
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The  results  of  a  previous  study  of  the  concept  of  a  stirred  reactor 
laser  with  lasing  restricted  to  the  zeroth  fundamental  band  are  compared 
with  the  predicted  power  output  of  a  multilevel  laser  for  the  hydrogen- 
fluorine-helium  (H2~F2-He)  system.  Lasing  is  assumed  to  occur  on 
all  fundamental  bands  for  which  the  calculated  maximum  gain  equals 
the  cavity  threshold  value.  Power  output  is  significantly  increased 
over  the  zeroth  band  calculations,  and  the  effect  of  cavity  pressure  is 
shown  to  be  less  pronounced  because  of  the  cascade  enhancement  from 
higher  bands.  Vibrational  de-excitation  caused  by  collisional  processes 
appears  to  limit  this  kind  of  chemical  system  to  low  cavity  pressures. 
The  effect  of  cavity  losses  was  evaluated  with  three  assumed  values 
for  the  cavity  loss  coefficient.  These  values  were  estimated  as  a 
small  fraction  (0.1)  of  the  average  small  signal  gain  coefficient  at 
different  cavity  residence  times. 
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